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PREFACE 

The  results  given  in  this  report  are,  in  large  part, 
identical  to  those  presented  in  number  IV  of  the  original 
series  of  mimeographed  lecture  notes  on  magneto-hydrodynamics 
issued  in  195i+»   The  class  of  parallel  flows  considered  in 
Section  2  of  this  report  has  been  extended  and  Section  5 
has  been  added. 
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Fluid  Dynainlcal  Analogies 

A»  A«  Blank  and  Harold  Grad 

1,   Introduction 

Although  the  combined  fluid  and  electromagnetic  equations 
generally  represent  a  far  more  complex  system  than  that  of 
classical  fluid  dynamics,  there  are  a  number  of  important 
instances  for  which  the  two  systems  are  mathematically  identical 
upon  appropriate  identification  of  symbolSo   The  existence  of 
these  exact  mathematical  analogies  permits  the  immediate  applica- 
tion of  a  wealth  of  known  resxilts  in  fluid  dynamics  to  the 
relatively  new  and  unexplored  field  of  fluid  magnetics.   One 
of  these  analogies  is  sufficiently  ramified  to  be  treated 
separately.,    The  results  given  in  this  report  are  misleadingly 
brief;  they  allow  the  application  of  many  known  techniques  and 
results  to  appropriate  problems  in  fluid  magnetics. 

The  fluid  dynamical  aspect  of  one  of  these  analogies,  i.e.  two- 
dimensional  non-steady  flow,  includes  the  major  portion  of  the 
entire  literature  of  fluid  dynamics.   On  the  other  hand,  the 
mathematically  identical  fluid  magnetic  aspect  is  much  more  of  a 
special  case.   This  remark  illustrates  the  far  greater  complexity 
of  the  totality  of  fluid  magnetic  theory  as  compared  to  conventional 
fluid  dynamical  theoryo 
MH-IX,  Equilibrium  and  Force-Free  Fields. 
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is 


The  general  set  of  equations  which,  is  \ander  consideration 
.2 


(1.1)  ft  "^  P  ^^^  u  =  0 

(1.2)  P  Hi  "^  S^^^  P  =  -  cxirl  B  X  B 

A 

(1.3)  p  ll  +  p  div  u  =  0,  or  §1=0 

(1.1].)         II  -  curl  (u  X  B)  =  0,  div  B  =  0 
together  with  the  auxiliary  relations 

(1.5)  E  =  -u  X  B 

(1.6)  J  =  -  curl  B 

(1.7)  q  =    X   div  E, 

We  use  the  notation  d/dt  to  denote  the  Lagrangian  derivative 

The  sound  speed,  a,  is  defined  by 

(1.9)  ^^  =^  »     P  =  ^^P>V 

A 

where  p  is  a  given  function  of  p  and  of  *^  (the  entropy  per  unit 
mass).   For  an  incompressible  fluid,  we  replace  (1,1)  by 

(1.10)  |£  =  0,   div  u  =  0 
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See  MH-VI,  Fluid  Magnetic  Equations  -  General  Properties, 
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and  Ignore  (1.3).   The  fluid  is  homogeneous  incompressible  if 
p  =  constant  throughout  rather  than  constant  on  each  particle 

A 

path.   Similarly,  if  ^  =  constant  throughout  tne  fluid  we 
call  the  flow  isentropic,  vjhereas  unoer  the  condition  (1.3)  the 
flow  is  called  adiabatico 

In  these  equations  we  have  assumed  infinite  conductivity 
and  have  omitted  displacement  ctirrent  and  electrostatic  forces. 
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2,   Parallel  Flows 

The   moinentum  equation    (1.2)    can  be  written 

(2,1)  P  1^  +   p(u.V)u  +   grad(p  +   ^^/Z\i)    =  -(B- V)B, 


The  form  of  this  equation  suggests  that  we  look  for  steady  flows 
in  which  B  is  parallel  to  u, 

(2.2)  B  =  Au; 

in  addition,  we  consider  only  an  incompressible  flow,  div  u  =  0, 
from  which  it  follows  that  \    is  constant  on  each  streamline 
(magnetic  line) 

(2.3)  H  =  (u-V  )A  =  0.    ■   , 

Since  equation  (l«l4-)  is  satisfied  identically,  all  that  remains 
is  to  satisfy  equation  (2,1)  which  takes  the  form 

(2.1^)     p"(u'  V)u  +  grad  p'"  =  0 

where   p"  is  a  constant  on  each  streamline, 

(2.5)     p'"'  =  P  -  \^/[i 

and  p"  is  given  by 

We  recall  that  (2,1^.)  together  with  div  u  =  0,  dp"/dt  =  0 

"  7  - 


(2,6)     p"^''  =  p  +  6^21 


defines  an  incompressible  steady  flow  with  density  p''  and 
pressure  p  ,    To  every  such  flow,  there  exists  a  steady, 
incompressible,  parallel  magnetohydro dynamic  flow;  both  p  and 
X  can  be  arbitrarily  specified  on  each  streamline  in  the  magnetic 
analogue.   Bernoulli's  law  holds,  viz.  the  quantity 

(2.7)    p"  +  ^"u  =  p  +  ^  pu 

is  constant  on  each  streamline. 

The  stability  of  the  degenerate  case  p   =  0,  p   =  constant, 
Tinder  the  homogeneity  condition  p  =  constant  has  been  treated 
by  Chandrasekhar  (proc.  Nat.  Acad.  Scl.,  l^,    273  (1956)\ 
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It  should  be  remarked  that  no  equation  of  state  is  required 


between  p"  and  p". 
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B,  z 


3,   Two-Dlmensional  Transverse  Flows 

¥e  consider  flows  in  which  all  quantities  are  independent 

of  z  but  may  depend  on  x,  y,  and  t.   The  magnetic  field  is 

restricted  to  the  z-direction, 
and  the  velocity  is  two-dimen- 
sional in  the  plane  (x,y).   It 
is  convenient  to  use  two-dimen- 
sional vector  notation,  interpreting 
the  single  component  of  B  as  a 
scalar »   Equation  (Idi-)  then  takes 
the  form 


-r- 


-^  1 


/ 


/ 


u 


Pig.  1 


(3«1)    dB  +  g  ^^^  ^  ^  Q_ 
dt 


In  other  words,  conservation  of  fliix 
takes  the  same  form  as  conservation  of  mass  in  this  geometry. 
The  reason  is  intuitively  clear;  in  a  plane  flow,  the  element  of 
area,  dxdy,  across  which  the  flux  is  conserved  is  the  same  as  the 
two-dimensional  "volume"  element  on  which  mass  is  conservedo 
Equation  (3.1)  combined  with  (1,1),  yields  this  result  in  the 
form 

(3.2)    ^(B/p)  =  0; 

the  ratio  of  B  to  p  is  constant  following  an  element  of  fluid. 
The  equations  (1,2)  and  (lo3)  may  be  rewritten  in  terras  of  the 
"total  pressure" 


(3.3) 


p  +  B  /2ti. 


and  "total  energy" 


(3.1;) 


=  e  +  B  /2p,p; 


namely,  noting  that  cvlpI  B  /  B  =  (B.V)B 
(B«  V)B  =  0  in  this  geometry,  we  have 


V(|-  B^)  and  that 


du 
dt 


(3.5)     P  ^  +  grad  p"  =  0 
and 


(3.6)     p  II  +  p  div  u  =  0; 

(the  entropy  form,  d  ij /dt  =  0,  can  also  be  used  instead  of  (3.6)  )» 
Despite  the  fact  that  equations  (3.5)  and  (3.6)  are  identical 
to  the  ordinary  fluid  equations  with  p"  and  e"  taking  the  place 
of  p  and  e,  this  is  not  an  exact  mathematical  analogue  because 
of  the  absence  of  a  conventional  equation  of  state.   More  precisely, 

the  complete  system, 

r 

^  +  p  div  u  =  0 
dt    ^ 

p  1^  +  grad  p'"  =  0 


(3.7) 


^    =   0 
dt     ^ 


dT 
\^   dt 


=   0  , 


)f  =  B/p 


has  two  quantities,  Y|  and  n'  =  B/p  which  behave  like  entropy, 
and  the  "thermodynamics"  has  three  degrees  of  freedom  instead 
of  two;  e.g.  p"  is  a  known  function  of  p,  r)  ,  and  Y)"' ,      An  identity 
which  will  be  fo\ind  useful  below  is 
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(3.8)    Td>^  +  S  dyf'  =  de'''  +  p"d(^) 


The  analogy  can  be  made  exact  by  imposing  additional 
restrictions: 

(a)  Suppose  Y)      -   constant  throughouto   B  can  be  eliminated  from 
the  equation  of  state  (3e3)  to  obtain 

(3.9)  p'''(p,  ^i  )  =  p(p,  ^  )  +  Ijl  P^  (*(*)^ 

This  equation  of  state  satisfies  the  usual  thermodynamic 
convexity  requirements;  in  particular  we  have 

(3.10)  Tdy^   =  de'""  +  p''''d(i) 

and 

(3.11)  (a''-)2  =  |e:  =  a^  +  i  p(yf'^-)2  =  a^  +  B^/t^P  • 

(b)  Suppose  fl   =  constant  throughout »   We  interpret  Ti     as  the 
thermodynamic  entropy  and  obtain  the  equation  of  state 

(3.12)  p'''(p,  ^)   =  p(p)  +  ^  p2(vj^*)^. 
Prom  the  relation 

(3.13)  de'"  +  p'"d(i)  =  ^  d/", 
we  see  that,  by  taking 

(3.11;)   t''  =  B/p., 

we  obtain  a  conventional  thermodynamic  structure.   The 
sound  speed,  a",  is  given  by  the  same  formula,  (3.11),  as 
previously^ 
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(c)   Suppose  that  Y)     =  f()?)«   More  precisely,  suppose  that  the 
curves  Vl(x,y)  =  constant  coincide  with  the  curves  h'  (x,y)=constant 
at  the  initial  instant;  we  describe  this  by  saying  that  m"  is  a 
(possibly  multi- valued)  function  of  V    initially*   Since  the  values 
of  V)    and  |o*  are  carried  with  the  fluid  element,  the  relation 
■)^  =  f(>7)  persists  for  all  time*   Either  n    or  Yj    can  be  chosen 
as  the  entropy;  let  us  take  y>  .   The  pressure  equation  of  state  is 

(3.15)  v''ip>Yi}   =p(p,  »j)  +^p2(f(>^))2  , 

the  sound  speed  is  as  before,  and  the  thermodynamic  identity 

(3.16)  T'"dn  =  de'"'  +  p'""d(i) 
follows  when  we  take 


(3.17)    T^"  =  T  +  i  p  f(»^)f'(y|) 


(d)   In  a  one-dimensional  flow  (variables  x  and  t),  the  relation 
rf    =   f{Y)  )  is  satisfied  for  arbitrary  initial  data,  and  the  formu- 
lation (c)  follows. 

All  conventional  two -dimensional  gas  dynajnical  theory  applies 
to  the  above  magnetic  analogues,  provided  that  no  specific  assTomp- 
tlons  have  been  made  regarding  the  equation  of  state.   One  of  the 
interesting  featiires  of  these  hydromagne tic  flows  would  seem  to 
be  the  possibility  of  creating  arbitrary  equations  of  state  at 
will  by  the  proper  choice  of  initial  conf iguratlon.   Another  way 
of  looking  at  this  is  that  one  can  construct  a  fluid  with  a  sound 
speed  which  is  an  arbitrary  function  of  position  much  more  easily 
than  in  a  conventional  gas, 
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[j..   Conservation  of  Circulation  and  Bernoulli's  Law 

Since  the  previously  described  analogies  are  exact,  the 
conventional  results  on  integrals  of  the  motion  follow  \ander  the 
same  conditions  as  quoted  for  the  conventional  fluldso   However, 
it  is  illuminating  to  derive  these  results  directly,, 

The  circulation  assigned  to  a  closed  curve  P  is  defined  by 

(1^.1)    C  =  d)  u.dx  =  qp  uo  ^  dcr 


r  r 


vjhere  cr  is  a  parametrizatlon  of  Po   For  a  ciirve,  P(t),  which 
is  moving  with  the  fluid  {cr     is  fixed  to  a  fluid  element)  we  have 


dC  _  /h  du  dx  ,^  ,  A  „„  du  ,_ 
TT  —  Qj  TT  '  J—  d(T  +  Qi'  u<*  -T—  dcT 
dt    T  dt   do-     1    dcr 


Y 


du 


^       1  2  ^    /  du 


(l;.2)       =  <©^  •  dx  +  0  dC^u")  =  (D^  •  dx 


dt 


=  (^|-|~(B»V)B  ^Vp'^'j  "    dx 


This  will  not  be  zero  and  yet  admit  any  significant  generality^ 
except  In  the  two-dimensional  transverse  flows  Just  treated,, 
In  that  case  (B*  V)B  =  0  and 


(^.3)    §     =     '^      ^^P^ 

P(t) 


Circulation  will  be  conserved  if  p  =  constant  (homogeneous 
incompressible  flow)  or  if  p'"'  is  a  function  of  p  alone.   The 
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latter  will  be  the  case  when  both  )?  and  h"  are  constant  through- 
out the  flow  (cf.  equations  (3.9),  (3.12),  and  (3.15)).   The 
general  condition  for  conservation  of  circulation  (cf,  equation 
(3.8))  is  that  a  potential,  _n_>  exist  satisfying 

ik.k)  Tdvj   +  ^  di^"   =  d(e'"  +  p'Vp)  -  -idp"  =  d_Q_ 

If  the  flow  is  irrotational, 
(1^.5)     u  =  V^ 


and  if  the  potential  _(^2_  exists  (e.g.,  if  )7  and  Y)"    are  constant), 

1 
2' 


1  2 
then,  using  (u.  V)u  =   curl  uXu  +  V(-5-u  ),  we  find  that 


or  (Bernoulli's  law) 

where  the  "total  enthalpy"  h'  is 
(lj..8)    h"""   =   e""'  +  p""'/P' 

For  a  steady  irrotational  flow  with  a  potential  _0_,  we  have 

(l|.9)     Vi^n)    +  ^Vp'''  =  0, 

from  which,  taking  the  dot  product  with  u  and  using  -rrr  =  w  V, 

we  obtain 

(I+.IO)     ^i^^   +  h'-''  -  XI)  =  0. 

In  this  case,  ^u  +  h"  -  _("}_  ^^  constant  on  each  streamline. 
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5,   One-Dimensional  Transverse  Flows 

This  geometry  is,  in  a  sense,  the  dual  of  the  one  treated 
in  Section  3.   We  assume  that  all  quantities  are  independent  of 

X  and  y  (the  relevant  variables  are  z 
and  t),  that  B  has  two  components,  B 
and  B  ,  and  u  has  only  a  z-component. 
We  again  have  (B-  V)B  =  0,  and  can  write 
^  the  momentum  equation  in  the  form 

^'   at    oz   p  oz 

The  equation  for  B,  interpreting  u  as 


For  the  magnitude  of  B  we  obtain 


(5.3) 


at 


1B|  +  ^(ulBl)  =  0. 


This  together  with  (5.1)  and  the  conservation  of  mass  is  a 
determined  system  (only  the  magnitude  of  B  enters  in  p").   More- 
over, it  is  exactly  the  same  system  that  was  treated  in  Section  3, 
specialized  to  one  space  dimension,  z,  instead  of  two,  (x,y)c   To 
complete  the  solution  of  the  original  problem,  we  note  that  (5.2) 
implies  that  the  direction  of  the  vector  B  is  a  constant  follow- 
ing a  fluid  element. 
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